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Abstract 

The von Karman-Howarth equations are derived for three-dimensional 
(3D) Hall magnetohydro dynamics (MHD) in the case of an homoge- 
neous and isotropic turbulence. From these equations, we derive exact 
scaling laws for the third-order correlation tensors. We show how these 
relations are compatible with previous heuristic and numerical results. 
These multi-scale laws provide a relevant tool to investigate the non- 
linear nature of the high frequency magnetic field fluctuations in the 
solar wind or, more generally, in any plasma where the Hall effect is 
important. 

Turbulence remains one of the last great unsolved problem in classical 
physics which has evaded physical understanding and systematic description 
for many decades. For that reason, any exact results appear almost as a 
miracle. In his third 1941 turbulence paper, Kolmogorov found that an 
exact and nontrivial relation may be derived from Navier-Stokes equations - 
which can be seen as the archetype equations for describing turbulence - for 
the third-order longitudinal structure function (Kolmogorov, 1941). Because 
of the rarity of such results, the Kolmogorov's four-fifths law is considered 
as one of the most important results in turbulence (Frisch, 1995). 

The derivation of the Kolmogorov's law uses earlier exact results found 
by von Karman and Howarth in 1938 (von Karman and Howarth, 1938): it 
is the well-known von Karman-Howarth (vKH) equation that describes the 
dynamical evolution of the second-order correlation tensors. Very few exten- 
sions of such results (vKH equations and four-fifths law) to other fluids have 
been made; it concerns scalar passively advected (Yaglom, 1949), such as the 
temperature or a pollutant in the atmosphere, and astrophysical magnetized 
fluid described in the framework of MHD (Chandrasekhar, 1950; Politano 
and Pouquet, 1998; Politano, Gomez and Pouquet, 2003). The addition in 
the analysis of the magnetic field and its coupling with the velocity field ren- 
ders the problem more difficult and, in practice, we are dealing with a couple 
of equations. 
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Signatures of turbulence in astrophysical flows are found in the solar 
wind (Goldstein and Robert, 1999; Matthaeus et al., 2005), the interstellar 
(Elmegreen and Scalo, 2004ab), galactic and even intergalactic media (Gov- 
oni et al., 2006). In the case of the interplanetary medium, we have access 
to very precise in situ measurements which show, in particular, the existence 
of a steepening of the magnetic field fluctuations spectrum at frequencies 
higher than 1Hz (Bale et al., 2005; Smith et al., 2006) whose origin may be 
attributed to nonlinear Hall-MHD processes (Galtier, 2006ab; Galtier and 
Buchlin, 2007). Efforts from observers are currently made to show the pres- 
ence of intermittency at high frequency. In this quest, any theoretical or 
model predictions about moderate or high-order correlation tensors is par- 
ticularly important for the understanding of solar wind, and more generally 
speaking, astrophysical turbulence. 

In this Letter, we derive the vKH equations for 3D Hall-MHD fluids. From 
these exact results, we give the equivalent of the Kolmogorov's four-fifths law 
for the velocity, magnetic and current density field correlations. 

We start our analysis with the following 3D incompressible Hall MHD 
equations 

(d t + v-V)v = -VP* + b ■ V b + uAv , (1) 
(d t + v- V)b = b • V v - d/V x [( V x b) x b] + r/Ab, (2) 

with V • v = 0, V • b = 0. The magnetic field b is normalized to a velocity 
(b — > ^/fi nmi b, with the ion mass and n the electron density), v is the 
plasma flow velocity, P* is the total (magnetic plus kinetic) pressure, v is 
the viscosity, rj is the magnetic diffusivity and dj is the ion inertial length 
(dj = c/ujpi, where c is the speed of light and ui pi is the ion plasma frequency). 
Equations ([I])-© may be rewritten more compactly as 

d t Vi = -<9jP* + bidfbi - vidfVi + vd^Vi, (3) 
dA = b e d e vi - Vfdfbi + d^Jgdibi - bid e Ji) + rjd^bi, (4) 

where J = V x b is the normalized current density. Note the use of the 
Einstein's notation. We see immediately that the third-order tensors that will 
appear in our analysis will be a combination of the velocity, the magnetic field 
and the current density. This makes an important difference with Navier- 
Stokes fluids for which the tensor used to derive the 4/5- law is built with the 
same (velocity) field. As shown below, it will have a direct impact on the 
kinematics. 
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Before deriving the vHK relations for Hall-MHD, one needs to introduce 
the kinematics adapted to this problem. The second-order correlation ten- 
sors, in the full isotropic and homogeneous case, may be written as (Batch- 
elor, 1953) 

R x (r) = (X i (x)X j (*!)) = F x r t r j + G x 6 ij , (5) 

where x' = x + r and X = v or b. F x and G x are four arbitrary functions of 
r 2 which will be specified later. The divergence free condition, d rj R x (r) = 0, 
on the velocity and the magnetic field leads to the relations 

4F X + rd r F x + r^a.G 1 = , (6) 

which will be used later. We now introduce the longitudinal and lateral 
functions as, respectively, 

R x = X 2 f x (r) and Rf ± = X 2 g x (r) . (7) 

The reference direction is the vector separation r such that, for example, the 
parallel component is the one along r. The correlation functions f v,b and g v,b 
are mainly decreasing (see [13] for more details about the function g) and 
satisfy the condition f v ' b (0) = g v ' b (0) = 1. From relations ©-([7]), we obtain 

itS(r) = X 2 (/% + C-5 t3 - r ~§)d r f x ) . (8) 

The third-order correlation tensors that will appear in our derivation are 
(Batchelor, 1953) 





J) 


= (vi(x)vj(x)v k (y^)) 

= AinrjVk + B x (riSjk + rj5 ik ) + £>ir fc <% , 


(9) 






= (^(x)6i(x)u fe (x')) 

= A 2 r i r j r k + B 2 (riS jk + rjS ik ) + D 2 r k S ij , 


(10) 






= («<(x)6i(x)6 fc (x / )) 

= A^iTjrk + B 3 ri5 jk + C 3 rj5 ik + D 3 r k 5ij , 


(11) 






= (J i (x)6 j (x)6 fc (x')) 

= A 4 rirjr k + B 4 ri5 jk + C A Tjb ik + D 4 r k 5ij , 


(12) 



where A m , B m , C m and D m are arbitrary functions of r 2 . Note that the 
last two tensors are not symmetric in the suffixes i and j which makes an 



3 



important difference with Navier-Stokes ffuids where only the velocity field 
is used to build the third-order correlation tensor. The direct consequence 
is that we need not three but four arbitrary functions to define initially 
these tensors. In the same way as before, we use the continuity condition, 
d rk S^ k (r) = 0, to constrain our system (Batchelor, 1953); it leads to the 
following relations 



rd r A 1)2 + 5A 1>2 + 2r- 1 d r B h2 = 


o. 


(13) 


rd r Di, 2 + 3D li2 + 2B 1>2 = 


o, 


(14) 


rd r A 3A + 5A 3A + r~ l d r B 3A + r~ l d r C 3A = 


o, 


(15) 


rd r D 3A + 3^3,4 + B 3A + C 3A = 


0. 


(16) 



Additionally, we note that S™ h (r) = whatever the value of m is (since it is 
a solenoidal first-order isotropic tensor); it gives the relations 

A lj2 r 2 + 25 1>2 + 3D 1>2 
A 3A r 2 + B 3A + C 3A + 3D 3A 



0. 
0. 



(17) 
(18) 



We introduce now the basic functions which include the parallel and perpen 
dicular components of the fields. We have 

y1,2 



fl,2 



-L-MI 
'1,2 
1 1 _l l_ 

3,4 



5 |,„ 
3,4 



_l l_ 1 1 

3,4 
±± 
3,4 



A h2 r 3 + (2S lj2 + D ia )r = Y 1>2 K li2 {r) , 
D li2 r = Y lj2 h 1:2 (r) , 
Bi, 2 r = Y li2 q 1>2 (r) , 

A 3A r 3 + (B 3A + C 3A + D 3A )r = Y 3A K 3A (r) 
D 3A r = Y 3A h 3A (r) , 
B 3A r = Y 3A q 3A (r) , 
C 3A r = Y 3A s 3A (r) , 



(19) 
(20) 
(21) 
(22) 
(23) 
(24) 
(25) 



where K m , h m , q m and s m are odd scalar functions and Y\ — v , Y 2 — Y 3 — 
vb 2 , Y4 = Jb 2 . Conditions f|T3|) - f|T6|) and f|T7|) - f|T8l simplify the expression of 
the third-order tensors, which write finally as 



+ 



Yl,2 

'2K li2 + rd r K li2 
4r 



' K h2 -rd r K^ 
2r 3 

{riS jk + rj5 ik ) 



riTjTk 
K 



(26) 



'1,2 



2r 



4 



Sijk( r 



Y, 



3,1 



K$A - rd r K 3A 



+ 



2r 3 

^3,4 + ^r^3,4/2 - q 3 ,4 



, 93,4 , 



(27) 



The first main goal of this Letter is the derivation of the vKH equations 
for 3D Hall MHD fluids. We start with equations (J3j) (SI) and use previous 
relations derived from the kinematics to find 

d t Rl{r) = {v^) + tydtVi) (28) 
- (v'jVtdeVi) - (v'jdiP) + vfa&^v'j) + vtydfai) > 



^(r) = (bidtfy + (b'jdtbi) (29) 
= (hb 1 ^) - Mdft) + (VjbidtUi) - (l/jVAbi) 
+ djdbiJ'Abj) - (bib'^Jj) + (bjJedth) 

After simple manipulations where we use, in particular, the divergence free 
condition and the homogeneity assumption, we get 

dtRij = d re [s]^ + Sj u — S?tj — Sj H ^j + 2vdr m Rij, (30) 
dtRij — d ri (^Sfj i — Sj U + S^j — Sf £ j + djSj ei (31) 
— diS^ji + djSffj — djSftj^J + 2r]dr ere Rij ■ 

Note that the pressure terms are suppressed because of isotropy (Batchelor, 
1953). These general dynamical equations reduce to a simple form for the 
diagonal part of the energy tensor, 

®tRii = 2d re - 5f fi ) + 2udr en Rii , (32) 
9tR b u = 2d re (Sea - Sf ei + di(Sf ei - S 1 ^)) + 2r]d^ ere R b u . (33) 

It is the basic equations from which it will be possible to derive the equivalent 
of the vKH relations. The introduction of $25]) and (J27J) into ([32]) and (JMD 
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gives 



(f 3 vb 2 \ 

—r e (A + rd^K, - — r £ (4 + rd r )K 2 J 



+ 2^<, 



2 r>v 

Avb 2 



d t R% = d r A — r e (K 3 + rd r K 3 /2-2q 3 ) 



+ ^^r e (K A + r0 r tf 4 /2 - 2q 4 )) + 2^ r ,^ 



By introducing 

K n + r9 r K n /2 - 2g ri 



K = d (r K ) K 

r 4 r 

for m = (1, 2) and n = (3, 4), we obtain 

3t((3 + r0 r )/V) = v^nK^-vtfd^nKz) 

+ 2vd 2 rere ((3 + rd r )rv 2 ) , 

d t ((3 + r0 r )/ 6 6 2 ) = 4rf 7 J6 2 ^(r^ 4 ) - 4^ 2 <9 r ,(r^ 3 ) 

+ 2^((3 + r9 r )A 2 ) . 

By noting the following identity (for isotropic turbulence) d 2 m 
(2/r)d r , and the more subtle relation 

(d 2 rr + \)(3 + rd r ) = (3 + rd r )^d r (r%) , 

we finally obtain after some simple manipulations 

d t ((3 + rd r )f v v 2 ) = v 3 (3 + rdr)^ 

-vb 2 {3 + rd r )K 2 + 2v (3 + rd r ) \d r (r A d r f v v 2 ) , 

r 4 

d t ((3 + rd r )f b b 2 ) = UiJb 2 (3 + rd r ) K A 
-Avb 2 (3 + rd r ) K 3 + 2r](3 + rd r ) \d r (r 4 d r f b b 2 ) . 
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A first integral of these equations is 

d t f v v 2 = v 3 K 1 -vb 2 K 2 + ^d r {r%f v v 2 ), (42) 

d t f b b 2 = -Avb 2 K 3 + ^d I Jb 2 K 4 + ^ld r (r 4 d r f b b 2 ). (43) 

These exact equations are the vKH relations for 3D Hall MHD. It is the first 
main result of this paper. 

We shall derive now the equivalent of the four-fifths law found by Kol- 
mogorov for Navier-Stokes fluids (Kolmogorov, 1941). We note the relation 
Rf\\{ r ) — (X 2 ) — !-Sfjj( r )> w ^ n general form for the structure function 

B$(r) = ((X^') - X^XjW - XjW)} , (44) 

where X = (v, b). Introducing this relation into the vKH equations (j4"2l - (j4"3~|) . 

we get 

d t {vf)-±d t Bf l{l = 
+ 

d t {b\) -\d t B\^ = 
+ 

and eventually 

B v 

dt(v\)-d t ^ = v^-vb^-^dr^drB^), (47) 
B b 

dtib^-dt^ = ±d I Jb 2 K±-\vb 2 K z -^d T {r^d r B\{). (48) 

We define the mean (total) energy dissipation rate per unit mass, e T , for 
isotropic turbulence, as 

$(«} + &}) = -(2/3)e T . (49) 

Exact scaling laws for third-order correlation tensors may be derived from the 
previous relations (I4"7T) - (|4"8"1) by assuming the following assumptions specific 



v 3 K x - vb 2 K 2 
2 ^d r (r%{{v 2 ) - l -B^) , 
Ad^b 2 ^ - Avb 2 K 3 
2 Jd r (r%.((bl) - Is},,)) , 



(45) 



(46) 
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to fully developed turbulence (Frisch, 1995). We first consider the long time 
limit for which a stationary state is reached with a finite e T . Second, we take 
the infinite (magnetic) Reynolds number limit (1/ — > and 77 — > 0) for which 
the mean energy dissipation rate per unit mass tends to a finite positive limit. 
Therefore, in the inertial range, we obtain at first order the relations 

- -e T r = v 3 (K x + rd r KJA) - vb 2 (K 2 + rd r K 2 /A) 
6 

- vb 2 (K 3 + rd r K 3 /2-2q 3 ) (50) 

+ d T Jb 2 (K 4 + rd r K 4 /2 - 2q 4 ) , 

which can also be written as 

~r Tr = ( s \\±± + \ s \\ II ll) ~(S\\±± + lS\\ ll 11) ( 51 ) 

The last step consists in introducing structure functions which gives, after 
some manipulations, the final result 

-Vr = B^+B^-2B^-Ad I (Sl ii -St lli ) 7 (52) 
= B^ + B$-2B$ + 4d I ([(Jxb)xby, 

with B$ = ((o{ - - ftM - 7fe )>- 

Equations fH2l - fj43|) and fl52l) are the main results of this paper. The 
former equations are exact for homogeneous and isotropic turbulence and 
the latter assumed additionally the existence of a large inertial range on 
which the total energy flux is finite and constant. The most remarkable 
aspect of these laws is that they do not only provide a linear scaling for 
the third-order correlation tensors within the inertial range of length scales, 
but they also fix the value of the numerical factor appearing in front of the 
scaling relations. Another important remark is about the fields used to build 
the third-order correlation tensors. Indeed, the convenient variables are not 
only the velocity and magnetic field components but also the current density 
components. Note that attempts to find a simple expression in terms of only 
structure functions failed, and therefore relations fl52|) seem to be the most 
appropriate. A similar situation was found for example in MHD flows when 
the magnetic helicity is included in the analysis (Politano and Pouquet, 1998; 
Politano, Gomez and Pouquet, 2003). 
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The vKH equations (j42p — (j43p derived here in the framework of Hall MHD 
are compatible with the one derived by Chandrasekhar (1950) for MHD when 
the large-scale limit (dj — > 0) is taken. Note that some minor manipulations 
have to be made in (Chandrasekhar, 1950) to prove the compatibility since 
the notation are not the same (for example, we have P = vb 2 K^). As ex- 
plained above, the notation used here seems to be more suitable for Hall MHD 
which therefore requires for a better understanding a complete re-derivation 
of the dynamical equations. In the same way, when the large-scale limit is 
taken, relation (1521 is compatible with previous works (Politano and Pou- 
quet, 1998; Politano, Gomez and Pouquet, 2003), which are also compatible 
with Navier-Stokes fluids when additionally the magnetic field is taken equal 
to zero. 

The exact results found here provide a better theoretical understanding 
of Hall MHD flows. They show that the scaling relation does not change its 
power dependence in the separation r at small-scales if the statistical correla- 
tion tensor used is modified. The interesting point to note is the compatibility 
with previous heuristic and numerical results (Biskamp et al., 1996). Indeed, 
a simple dimensional analysis gives the relations r ~ b 3 for large-scales, and 
r 2 ~ b 3 for small-scales (since J ~ b/r), which give respectively the magnetic 
energy spectrum E ~ k~ 5 / 3 and E ~ k^ 1 ! 3 . Therefore and contrary to the 
appearance, the exact results found may provide a double scaling relation. 

These multi-scale laws provide a relevant tool to investigate the non- 
linear nature of the high frequency magnetic field fluctuations in the solar 
wind whose (dissipative vs dispersive) origin is still controversial (Goldstein 
et al., 1994; Markovskii et al, 2006; Stawicki et al., 2001; Galtier, 2006ab; 
Galtier and Buchlin, 2007). The use of multi-point data may give informa- 
tion about both the magnetic field and the current density which can be used 
to check the theoretical scaling relations. The observation of such a scaling 
law would be an additional evidence for the presence of a dispersive inertial 
range and therefore for the turbulent nature of the high frequency magnetic 
field fluctuations. The recent observation of the Yaglom MHD scaling law 
(Sorriso-valvo et al., 2007) at low frequency provides a direct evidence for 
the presence of an inertial energy cascade in the solar wind. The theoretical 
results given here allows now to extend this type of analysis to the high fre- 
quency magnetic field fluctuations and, more generally speaking, to better 
understand the role of the Hall effect in astrophysics, like e.g. for the magne- 
torotational instability in cool proto-stellar disks (Wardle, 1999; Balbus and 
Terquem, 2001), or in laboratory fusion plasmas. In such situations, isotropy 
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is often broken because of the presence of a strong large-scale magnetic field 
(see e.g. Muller et al., 2003) and a generalization of the present description 
to anisotropic turbulence is then necessary. 
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I would like to thank H. Politano for useful discussion. 
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